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Abstract 

We investigate superconducting systems with the use of the phenomeno- 
logical Landau's theory of second order phase transitions, including into 
the considerations the critical behaviour of the chemical potential. We 
derive in this way a variety of new thermodynamical relations at the 
critical point. Twelve basic relations connect critical jumps of different 
thermodynamical quantites (specific heat, chemical potential derivatives 
with respect to temperature, pressure (volume) and number of particles, 
volume (pressure) derivatives with respect to temperature and pressure 
(volume)) with the critical temperature or its derivatives with respect to 
the number of particles or pressure (volume). These relations allow to 
find plenty of cross-relations between different quantities at the critical 
point. The derived formulae can practically be used in many cases to find 
such thermodynamical quantities at the critical point which are extremely 
difficult to measure under the assumption that the other ones are already 
known. We additionally perform a test of the two derived relations by 
using two-band microscopic model, describing superconducting systems. 
We calculate the specific heat, order parameter and chemical potential as 
functions of temperature to show that the tested relations are very well 
fulfilled. 

Subject classification: 05.70.Fh; 74.20. De; 74.25. Bt 

Introduction First observations of the chemical potential critical behaviour 
at second order phase transitions has been performed for superconductors (see 
Refs [1-10]). It has been found that at the critical temperature Ts (phase tran- 
sition: superconductor - normal system) the temperature dependence of the 
chemical potential exhibits a kink (change in slope) , manifesting in this way the 
appearance of the phase transition (cf also Refs [11-23]). For superconducting 
systems (see Ref. [6]), an interesting relation between the jump of the chemical 
potential temperature derivative and the jump of the specific heat at the critical 
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point has been derived with the use of the phenomenological, Landau's theory 
of second order phase transitions (cf Refs [24, 25]). It is, however, possible to 
include in a more complete way the chemical potential critical behaviour into 
the general scheme of the Landau's theory for superconducting systems to derive 
new relations between different thermodynamical quantities at the critical point. 
Starting from the Gibbs potential we derive six basic relations and further six 
ones are resulting from the Helmholtz potential. There are two relations con- 
necting the critical jump of the specific heat at constant pressure (volume) and 
the critical temperature. The next two relate the critical jump of the chemical 
potential temperature derivative and the critical temperature derivative with 
respect to the number of particles (constant pressure (volume)). Critical jumps 
of the volume (pressure) derivatives with respect to temperature and pressure 
(volume) can be connected with the critical temperature derivatives with respect 
to pressure (volume), resulting in four relations. The next four relations con- 
nect critical jumps of the chemical potential derivatives with respect to pressure 
(volume) and the number of particles with the critical temperature derivatives 
with respect to the number of particles and pressure (volume). These twelve 
basic relations can be combined together by forming ratios of different critical 
jumps. In this way we eliminate an unknown parameter of the Landau's theory 
(second Landau's parameter) and we connect directly the mentioned ratios with 
the critical temperature of the system and its derivatives. These relations can 
be very useful in the practice because they allow to find an unknown critical 
jump (or unknown critical temperature derivative) when the other quantities 
are already known from experimental data. One can additionally eliminate the 
critical temperature derivatives from these relations to obtain six independent 
cross-relations between different critical jumps and the critical temperature of 
the system. These formulae can also be practically used to find an uknown 
critical jump under the assumption that the other critical jumps are already 
known. We additionally perform a test of two basic relations, desribed above, 
when independently calculate the second Landau's parameter from these two 
relations. To this purpose we applied two-band, microscopic model, describing 
superconducting systems. The numerical calculation of the specific heat jump 
and the jump of the order parameter at the critical point, as well as, the critical 
temperature derivative with respect to the number of particles (electrons) re- 
veals almost the same value for the second Landau's parameter when calculate 
it from these two independent relations. It seems to be very strange that the 
analytical derivation of the formulae, we have numerically tested, is extremely 
difficult to do for this particular, microscopic model (we could not perform it) 
because of its great complexity (to many implicit equations to solve) whereas the 
phenomenological derivation is not only general (independent from the assumed 
microscopic model) and straightforward but also true (!). 

Phenomenological relations Within the phenomenological theory of second 
order phase transitions (see Refs [24, 25]) one assumes the existence of the order 
parameter x. For simplicity, we assume that the system can be described by 
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only one order parameter |^. For superconducting systems the order parameter 
X = I Asp where is the superconducting order parameter (or its amphtude 
in the case of anisotropic superconductivity). The order parameter should fulfil 
a selfconsistent equation of the general form 

x^f{T,p,N,x) (1) 

where the variables T, p, N denote the absolute temperature, pressure and 
number of particles of the system. This equation can formally be resolved to 
give an equation 

x = xiT,p,N) (2) 

of such a property that a; = for T^Tg where Tg = Ts{p, N) is the critical tem- 
perature of the transition from superconducting to normal phase. One assumes 
also that the thermodynamic properties of the system can be desribed by the 
Gibbs potential 

Z = ZiT,p,N,x). (3) 

In the vicinity of T = (T < Tg) the order parameter x is very small and 
therefore we can expand the Gibbs potential (3) into the Taylor series (we 
retain quadratic terms) 

Z{T, p, N, x) - Zo [T, p, N) + a{T, p, N)x + ^b{T, p, N)x^ (4) 

where a and b are first and second derivatives of Z with respect to x at the 

dZ d^Z 

point a; = 0. The mininmm principle (-^ — 0, > 0), aplied to (4) results 

ox dx-^ 

in the expression 

x = -^ (6>0). (5) 

To obtain a general form of the phenomenological parameter a = a{T,p, N) in 
the vicinity of the critical point we can apply again the Taylor expansion with 
respect to a small quantity {T — Tg). We get 

aiT,p,N) = a{Tg + {T^Tg),p,N) 

« a{Tg,p,N) + {T-Tg)(^) . (6) 

V^^ / p,N,T=T^ 

However, at T = the order parameter x = Q.lt means also that a(Tg,p, N) = 
(see (6) and (5)). Thus, we obtain 

(T-Tg)(^) ,T^Tg 
a{T,p,N)={ / p.N,T=T^ (7) 

3 ,T>Tg 



^See also Note added in proof 
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and 

{T - T,) (da\ 



X = x{T,p,N) = { ^ \^^Jp,N,T=T, ' _ (8) 

,T>T, 
Eliminating the phenomenological parameter b from (4) by means of (5) we get 

Z = Zo + ^ax (9) 
and using the thermodynamic relations 

s=-m] . (10) 



dT 



So = -{'4^) (11) 



dT 



p,N 



we can write 



1 /9a\ 1 dx\ 



The specific heat difference, produced by (12) is given by relation 



p,N 

Tfd^\ _rp(da\ fdx_\ 



2 ydT"^) 



i2. . 

(13) 

p,N 



To obtain the jump of the specific heat Acp^N at the critical point we have to 
insert T = Tg into (13) taking also into accout that x = a = Oa,tT = Ts (see 
(7), (8)). Thus, we obtain 

A».» = -r,(^) m ^ (14, 

V'^-' / p,N,T=T, / p,N,T=T, 

However, when looking at (8) we can see that 

=-l(^\ (15) 



^ J p,N,T=Ts ^ / p,N,T=Ts J p,N,T=Ts 

fdx\ 

where A — — I is just the jump of the temperature derivative of the 

V^-' J p.N,T=T^ 

order parameter at T = T, because the order parameter x vanishes for T ^ Tg. 
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Thus, we can rewrite (14) with the use of (15) by ehminating ( 
We obtain 

r H ^ 

Ac,.N = hT, a(^) . (16) 

The chemical potential of the system can be introduced from the general 
expression for the Gibbs potential using the relations 



These relations applied to (9) result 

1 f da\ 



Mo 



dZo\ 



1 f dx 



(17) 



(18) 



Thus, the jump of the chemical potential temperature derivative at T = Tg 
(a; = 0, a = for T > Tg) is equal to 



A 



N,T=T, 



^ J p,N,T=Ts \ J p,N,T=Ts 

1 f da\ ( dx' 



1 f da 

2 \&f 



dx 

p,N,T=T, \^^^ / T=Ts,p 



Using (7), (8) and (15) we find (T^ = Tg{p,N)) : 
da\ _ _[dTs\ (da 



(19) 



p,N,T=T, 



and 



^ J p,N,T=T^ 



(20) 



dx \ 
dNj 



1 fdTA (da 



T=n,p b \dN J p\dTy 



Introducing (20) and (21) into (19) we get 



J p,N,T=T^ 



/ p,N,T=Ts 



-\ 2 



(21) 



(22) 
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It is also easy to find that 



.dp, 



T=Ts,N 



<^P / N 



da 
df 



p,N,T=T, 



(23) 



dp 



dp 



dT 



N \^-^ / p,N,T=Ts \ / N \^-^ / p,N,T=Ts 

(24) 

Wc can perform similar calculations as above to find critical jumps of the follow- 



ing thermodynamical quantities: A (f^)^ , A (^^^ , ^ 



Starting from the expression for the Gibbs potential (9) 
with the use of the thermodynamical relations /x = (fy)^,^ (Mo = {^t)tp^^ 
^ " (^)tjv " (^)tjv'' """^^ ^^"""^ formulae (15), (20), (21), (23) 



and A ( ^) 



and (24) we obtain 



d^ J p,N,T=Ts 



dp / N 



dT 



-\ 2 



p,N,T=T, 



(25) 



.9pJt=t„n \dNjp\dpyj^ 



A 



dT 



p,N,T=T, 



(26) 



and 



A 



dp ) T=Ts,N 



\dN 



{dp 



N. 



A 



A 



dx\ 

J p,N,T=Ts 

dx\ 

dT J p^N,T=Ts 

■ dS 



= -(^)..,,(m^_ = -{^] and 

p,N 



(27) 



(28) 



Using the MaxweU relations (§ )^_^ - var/p.jv' Vaw/r.p - \dT, 

{^)rpp — (^^^ we can see that the other critical jumps A 

, A and A (|^)j,^y can easily be obtained from (25), (22) and 

(26), respectively. In other words, there are only six independent expressions 
describing the critical jumps of the first derivatives ((16), (22) and (25)-(28)). 
These expressions are very important because with the use of them we can 
find many interesting cross-relations. Dividing e.g. (22) and (25)-(28) by (16) 
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wc obtain five independent ratios (the specific heat critical jump is finite for 
superconductors) : 



A ( 



N,T=T, 



Ac, 



'p,N 



T=Ts,N 



ON 



dp 



Ac 



'p,N 



A I 



T=T,,p 



Ac 



and 



'p,N 



Ac 



'p,N 



dlnT^ 
dN 

d\nT, 
dN 

dlnT^ 
dp 



N 



N 

dp 

dTs 
dN 

dls 
dp 



N 



(29) 
(30) 

(31) 

(32) 

(33) 



JV 



It is, however, possible to find also another ratios (equivalent to the above 

relations). When eliminate the derivatives and (^■j from the above 

expressions we can obtain three independent cross-relations connecting different 
critical jumps: 



9p J T=Ts,N 

A 



■ Ac, 



'p,N 



-T.Ai^ 

dT 



p,N,T=Ts 



.A(^) 

dT J 



dji 
dN 



^ 1 ■ Acp,jv = -T, 



T=Ts,p 



and 



\ 9p 



A(^ 

dT 



A(?S 



p,N,T=T, 
-\ 2 



p,N,T=T, 



dT 



-\ 2 



p,N,T=Ts 



(34) 
(35) 

(36) 



• Acp,jv = -Ts 

' T=T,,N 

These relations allow to find some critical jumps which are difficult to measure, 
knowing the other from experiments, without necessity to know the derivatives 
of the critical temperature. 

We can also derive several additional relations, when using the thermody- 
namical relationship resulting from the Gibbs potential (1) 



V = 



dZ^ 

dp 



V{T,p,N,x). 



(37) 



T,N 



The equation (37) can formally be resolved to obtain the equation of state 

p = p{T,V,N,x). (38) 
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We can insert it into (1) to obtain 



Z = Z{T,V,N,x) 



(39) 



which is nothing else but the Helmholtz potential. We can further apply exactly 
the same procedure as described above starting again with the Taylor expansion 
of the expression (39) with respect to the order parameter in the vicinity of the 
critical point. Therefore there is no need to demonstrate here all the details of 

the performed calculations. We restric;t ourselves to present only final results. 
The analogous formulae (cf (16), (22) and (25)-(28)) for the critical jumps, 
resulting from (39), have the form 



Ac 



V,N 



A 



A 



df 



dp 
df 



V,N,T=T, 



= -b 



A(^ 

dT 



OTA 
dNjy 

on 



-\ 2 



V,N,T=Ts 



A 



dT 



V,N,T=T^ 



V,N,T=Ts 



dV 



N 



\^-^ / V,N,T=Ts 



(40) 
(41) 

(42) 



A ^ =-b 



T=T,,N 



dN 



dTs 
dV 



N 



V,N,T=T, 



(43) 



and 



djj, 
dN 



= -b 



T=Ts,V 



dN), 



dx 
'df 



V,N,T=Ts 



(44) 



A^ 



dp 
dV 



T=T,,N 



dV 



N. 



A(^ 

dT 



V.N.T=Ts 



(45) 



V,N 



T,N 



The constant b has the same meaning as before but actually b = b{Ts ,V,N). To 
derive (40)- (45) we have applied thermodynamical relations S = 

V = - (# and , = (11)^^^ (^0 = - (#)^_^, .0 = 
and ^lo = (y^^^^ where Zq = Zo{T,V,N) in analogy to (4)). Because 
of the Maxwell relations (|f)^_^ = (t)^^^, (|^)^^^ = - (|^)^_^ and 
im)T,v = - ^"ti^^l j^^P« ^ iw)T=T.,N ' ^ and 



9Zo 
dV 



8 



A (^)^_^ ^ can easily be found from (42), (43) and (41), respectively, there- 
fore there are again only six independent formulae (40)- (45). When e.g. devide 
(41)-(45) by (40) we get: 



and 



V,N,T=T^ 



Acv,N V 9N 



V,N,T=Tg 



dlnT, 



N 



Acv,N V 9V 

T=T,,N f dlnTs\ f dT, 



Acv,N V dN Jv\dV 

^i^)T=n,v __ f dlnTA fdT^ 



N 



(46) 
(47) 
(48) 
(49) 

Acv,N V 9V j^KdVj^- 

There are, however, many other possilities to form the ratios, leading to equiv- 
alent relations. 

At last, when eliminate from (46)-(50) the critical temperature derivatives, 
we obtain (cf (34)- (36)) 



A 



dp 



dV 



Acv,N = TsA 



T=Ts,N 



df 



V,N,T=Ts 



dp 
df 



A 



Acv, 



N 



T=Ts,V 



df 



V,N,T=Ts 
2 



and 



A(^ 

dV 



Acv,N = Tg 



T=Ts ,N 



A(^ 

dT 



V,N,T=T, 



(51) 

(52) 

(53) 



Thus, the derived formulae (16), (22), (25)-(28) together with (40)-(45) al- 
lowed to find many useful thermodynamical cross-relations for superconducting 
systems at the critical point. 



9 



Numerical Test: superconductor - normal system transition In the 

following we perform a validity test of the formulae (40) and (41) when using 
microscopic model, describing superconducting properties of a system with two 
hybrydizcd electronic bands. The Cooper - pairs attraction, responsible for the 
superconducting properties, we choose in the form of the intersite interaction 
restricted to only one band (one superconducting order parameter). The Hamil- 
tonian of this model can be written in the form (great canonical ensemble) 

H = H-i^N (54) 

where 



a=l i,j,a i,a 

and 

^ = EE-S- (56) 

a=l i,a 

Two electronic bands are characterized by the hopping intergrals t^J (a = 1, 2; 

t^j^- = 0, t[^^ — t). The parameters V and Rij denote the hybridization and 
intersite Cooper - pairs attraction, respectively. The electronic annihilation 

(creation) operators (a = 1, 2; a = t, i) are denoted by c-"^-* (c-'^'' ), respectively, 
where i is the lattice site index. For simplicity, we consider simple cubic (sc) 
lattice. Similarly to Ref. [26] we apply the mean field approximation (see the 
third term in (55)) 

(2)+ (2)+ (2) (2) _ / (2)+ (2)+\ (2) (2) / (2) (2)\ (2)+ (2) + 

-(4'*4r)(cS«s)- 

Introducing (57) into (55) and applying Fourier transformation we obtain 

H = -E^k(cL^c^') + EE4"^(<?+-^ti) 

k Q=l k 

+E'^kcLl,cL^l+EA''4!lcLl,, (58) 

k k 
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^=EE(<T+"-ti) (59) 



a=l k 



where 



^k = ^E^(k-k')(cLt,cg), (60) 



k 



4^ = — ^^(k). (61) 

4^=i--^^^(k), (62) 
i?(k) = 2RoF{\i), (63) 

i^(k) = cos(fcxo) + cos(A;yo) + cos(A;za) (64) 

and Af is the number of lattice points. We assume that the bandwidth = W 
and W(i) = (5W ((5 < 1). In the following we assume that (cf Ref. [26]) 

Ak = A,F(k) (65) 

where As (wc assimic A* = Ag) is the amplitude of the superconducting order 
parameter (s* - symmetry). Looking at the expressions (61)-(63) and (65) we 
see that the only dispression in our system is generated by the formfactor (64) 

w 

2 

of the sc lattice. It means that we can convert jj- ^(...)^ J po{e){...)de where 

k _ 

2 

Po{e) is the true density of states of the sc lattice given by the Jelitto formula 
[27]. Applying the standard Green's function approach (cf e.g. [28]) we can find 
a set of transcendental equations of the type (a= 1,2; a = t,|) 

(n(")) = /(")(T,A2,^), (66) 
A2 = A^/(T,A^,^) (67) 

and 

2 

^^(n("))=n (68) 

q:=1 a 

(N) . 

where n = ^ is the average number of electrons per lattice site. The calcu- 
lations of all necessary quantities to check the formulae (40) and (41) are very 
laborious and lead to very complicated expressions (to long to present here). 
Therefore we are forced to describe only the general calculation scheme and 
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A*ltf[eV/K] 




(^T)„„*10P[eV/K] 



(C) 



=6*10"'[eV/K] 



87.0 87,5 88,0 88,5 89,0 87,0 87,5 88,0 88,5 89,0 87.0 87,5 33,0 88^5 89,0 

teiniieiaure [K] tmiierature [K] tanperalure [K] 

V10^[eV/Kl (d) PWIK^l (e) e^^/910,„-10=[(eV)^/Kl (f) 

OlnTj/9n),^1.81 0,0000- 



(Ac^,J,_,=32*10-' 
[eV/Kl 




87,0 87,5 SS,0 &&,5 89,0 0^ 0,9 1,0 1,1 1,2 15 87,0 87,5 33,0 83,5 39,0 

temperature [K] rt 

temperature [K| 



A(apy5T),^„_/(Ac,^„)„=i.8i 



Figure 1: Plot of the coefficients A (Fig. la), B (Fig. lb) and (|^) (Fig. Ic) entering 

V ^^^-^ J V,n 

into the specific heat per lattice site (Fig. Id) as functions of temperature in the vicinity 
of the critical point. The plots in Fig. le, f show the superconducting critical temperature 
derivative with respect to the average number of electrons per lattice site n as function of n 
and the temperature derivative of the order parameter as function of temperature. The 
model parameters: W= 3 eV, <5=0.25, Ro=0.6 eV, t=l eV, y=0.15 eV, n=0.97. 



to present final, numerical results. The implicit equations (66)-(68) allow to 
calculate the temperature dependence of the square amplitude (it is just 
the order parameter in the sense of the Landau's theory) and n as function of 
temperature. The calculation of the specific heat per lattice site cv,n can be 
performed with the use of (58) when calculating 

(^) = (69) 

and differentiating it with respect to T. This differentiation is, however, very 
laborious because and /x fulfil the equations (66)- (68) and are also temper- 
ature dependent. In other words, we have to differentiate them with respect to 
T to calculate the temeprature derivatives of A^ and fj, (necessary to check the 
formulae (40) and (41). The transition temperature Tg (phase transition from 
superconducting to normal state) can be obtained from the equation (67). It is 
easy to see that Tg should fulfil the following equation 

f{Ts,l^) = l. (70) 
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The function / in (70) is exactly the same as in the rhs of (67) but we 
have additionally assumed that the order parameter is equal to zero (T = T^). 

( dT \ 

We should also calculate the derivative J . It can be performed by the 

\On J y 

differentiation of the equations (70), (66) and (68) with respect to n (we put 
T = Ts therein (the order parameter is equal to zero)). In such a way we have 
all necessary quantities to check the formulae (40) and (41). These formulae 
we should first rewrite in the following way {N — > (N) = Afn, x = A^, 
b — * ^—jr (*^® second Landau's parameter per lattice site)). Thus, we obtain 



Acv.i 



T=T, 



= b 



(71) 



and 



V.n.T=T^ 



\ dn )^ 



)v,n. 



T=T, 



= b. 



(72) 



The specific heat per lattice site can be written in the form 



cv,n {T) = A{T) + B{T)[^ 



(73) 



V,n 



The results of our calculations are presented in Fig. 1. We see (Fig. la, c) 
that A (T) and the temperature derivative of the chemical potential exhibit 
jumps aX T = Tg whereas the coefficient B (T) is a continuous function at 
the critical point (Fig. lb). It results in the jump of the specific heat per 
lattice site cy_„ (Fig. Id) at T = T^. It would be interseting to know whether 
the temperature behaviour of the coefficients (Fig. la, b, c) in (73) is a general 
property of all superconductors (when write the specific heat in the form (73)) or 
only a property of this particular model. This question needs a verification from 
the experimental (nobody has proved it till now) and also from the theoretical 
point of view (when consider another microscopic models). From the other 
hand, the model we have applied here, should qualitatively reproduce all general 
properties of superconductors. Therefore it is very probably that only coefficient 
A and the temperature derivative of the chemical potential exibit jumps at 
the critical point whereas the coefficient i? is a continuous function for these 

dTs dA^ 

materials. The plots of the derivatives ( ^ ) and { ^ " ) are given in Fig. 

on V oT v,n 

le and Fig. If, respectively. Calculating the Ihs of (71) and (72) independently 
by means of the data displayed in the Figs Ic-f we obtain almost the same value 
of the second Landau's parameter per lattice site b « 8.6 eV~^. It also means 
that the formula (46) is automatically fulfilled (see Fig. Ic, d, e). 
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Conclusions The general, phenomenological relations (16), (22) and (25)- 

(28) together with their analogs (40)- (45) seem to be very important for super- 
conducting systems in the practice because they connect different macroscopic 
thcrmodynamical quantities at the critical point. They also allowed to obtain in- 
teresting cross-relations (29)-(36) and (46)-(53) between different critical jumps 
and the critical temperature derivatives or between the separate jumps. Due to 
a variety of relations wc have a possibility to estimate in an indirect way many 
of the critical jumps without performing experiments. Such a possibility can es- 
pecially be interesting in many cases where experiments are difficult to perform 
and the information gained in this way can be decisive as e.g. the indication for 
the hole- type superconductivity in the high-Tg superconductor Y Ba2Cu30r-s 
(see Refs [6], [7] for details). 

Note added in proof The derived critical relations (29)- (36) and (46)- (53) 
are also valid for a more general case of superconducting, correlated systems with 
many order parameters (multiband systems). As a consequence of electronic 
correlations the critical temperature Tg is the same for all order parameters (all 
of them vanish for T > Tg). To make a proof of the statement that the relations 

(29) -(36) and (46)-(53) are valid in this case we consider the Gibbs potential in 
the form 

n 1 ^ 

Z{T,p,N,xi,X2,...,Xn) = Zo{T,p,N) + Y,aiXi + ^ ^XiXj (N. 1) 

as a generalization of the Gibbs potential (4) after Taylor expansion with respect 
to the order parameters Xi (i = 1,2, ...,n). The parameters a^, bij = bji = 
1,2, ...,n) have similar meaning as in the formula (4). Using the orthogonal 
transformation 

n 

Xi = Y^ aijXj (N. 2) 

i=i 

the bilinear form in (N.l) can always be diagonalized {bij = bji) and instead of 
the formulae (N.l) we obtain the following expression (similar to (4)) 

n 1 ^ 

Z{T,p, N, xi,X2, Xn) = Zo{T,p, N) + Y, + 2 ^) 

i=l i=l 

where Oj = ^ and bi {i = 1,2, ...,n) are the eigenvalues of the matrix 

i=i 

bij. The quantities .t^ (z = 1, 2, n) can be treated as "new" order parameters 
because of the reciprocal relation 

n 

Xi = Yxjay, (N. 4) 
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and due to the fact (see (N.4)) that Xi {i = 1, ...,n) vanish at the same critical 
temperature as the order parameters Xi {i = l,...,n). From the minimum 
principle =0, i = 1, ...,n) we obtain 



and this equation can be used to eliminate bi {i = 1, ...,n) from the expression 
(N.3). We obtain in this way very similar expression to (9), it is 



To each and Xi {i — 1, n) we apply analogous assumptions to (7) and (8) 
where the lower index i should be added. Repeating the calculations similar 
to the described above (see (10)-(28)) we obtain analogous expressions to (16), 
(22) and (25)- (28). These new expressions, however, can easily be obtained 
from (16), (22) and (25)-(28) when replace b[A{^)p^N^T=Ts? on the rhs of 

n _ 

them by a new quantity ^i[^{^)p,N,T=Ts]'^ ■ This new quantity is irrelevant 



when forming the quotients (sec (29)-(33)) because it cancels out. Thus, the 
formulae (29)- (33) and the resulting cross relations (34)- (36) remain unchanged. 
Exactly the same can be said about the formulae (46)- (53) which remain also 
the same when apply the procedure described in this section with respect to the 
Helmholtz potential Z = Z{T, V, N, xi,X2, Xn) as a generalization of (39) for 
the case of superconducting systems with many order parameters. 
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